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Abstract 

This paper is a continuation of the work started in [1] and [2]. We extend one of the 
theorems that gave a way to construct equilateral triangles whose vertices have integer 
coordinates to the general situation. An approximate extrapolation formula for the 
sequence ET(n) of all equilateral triang les with vertices in {0, 1, 2, n} 3 (A 102698) 
is given and the asymptotic behavior of this sequence is analyzed. 

1 Introduction 

It turns out that equilateral triangles in Z 3 exist and there are unexpectedly many. Just 
to give an example, if we restrict our attention only to the cube {0, 1,2,..., 1105} 3 we have 
2,474,524,936,846,512 of them. In [1] it was shown the first part of the following theorem and 
the second part about the converse was only proven under the hypothesis that gcd(d, a) = 1 
or gcd(d, b) = 1 or gcd(d, c) = 1. The main result of this paper is to show that one can drop 
this condition. 

Theorem 1.1. Let a, b, c, d be odd positive integers such that a 2 + b 2 + c 2 = 3d 2 and 

gcd(a, b, c) = 1. Then the points P(u, v, w) and Q(x, y, z) whose coordinates given by 

u = m u m — n u n, [ x = m x m — n x n, 

m v m — n v n, and \y = m y m — n y n, (1) 
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with 



m x = — \[db(3r + s) + ac(r — s)]/q, n x = —(rac + dbs)/q 



m y = ^[da(3r + s) — bc(r — s)]/q, n y = (das — bcr)/q 
m z — (r — s)/2, n z = r 

and (2) 
m u = —(rac + dbs)/q, n u = — \[db(s — 3r) + ac(r + s)]/q 
m v = (das — rbc)/q, n v = \[da(s — 3r) — bc(r + s)]/q 
m w = r, n w = (r + s)/2 

where q = a 2 + b 2 and (r, s) zs a suitable solution of 2q = s 2 + 3r 2 which makes all the 
numbers in (2) integers, together with the origin (0(0,0,0)) forms an equilateral triangle in 
Z 3 contained in the plane 

V a ,b, c ■= {(a,f3,j)\aa + b(3 + 07 = 0} 



and having sides-lengths equal to d^2(m 2 — mn + n 2 ). 

Conversely, there exist a choice of the integers r and s such that given an arbitrary 
equilateral triangle contained in the plane V a ,b,c with one of the vertices the origin, the other 
two vertices are of the form (2) for some integer values m and n. 

The condition 

min(gcd(d,a), gcd(d,b), gcd(d,c)) > 1 (3) 
defines a, b, c, and d as a degenerate solution of the Diophantine equation 

a 2 + b 2 + c 2 = 3d 2 . (4) 

The first d that admits such degenerate decompositions in a, b, and c is d — 1105. There 
exactly seven of them: 

(a, b, c) G {(731, 1183, 1315), (475, 1309, 1313), (299, 493, 1825), 

(1027, 1139, 1145), (187, 415, 1859), (265, 533, 1819), (493, 1001, 1555)} 

Our study started with the intent of computing the sequence ET(n) of all equilateral 
triangles with vertices in {0, 1, 2, n}. These values were calculated by the first author 
for n < 1105 using a improved version of the code published in [2] and translated into 
Mathematica (see A102698 in The On-Line Encyclopedia of Integer Sequences). 

One of the parametrizations like in (2) in the case a = 731, b = 1183 and c = 1315 is 
shown below: 

P = (901m - 1428n, -1157m + 221n, 540m + 595n), 

Q = (-527m - 901n, -936m + 1157n, 1135m - 540n), m, n G Z. 
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For other values of r and s we get formally different parametrizations but they are all 
equivalent in the sense that they can be obtained from the above by changing the variables 
in the expected ways: 

(m, n) — > (— m, —n), (m, n) — > (m — n, m), (m, n) — > (n — m, n), 

(5) 

(m, n) — > (m — n, —n) and (m, n) — > (n — m, —m). 

These changes of variables leave invariant the quadratic form involved in the side-length 
of the triangle the various points given by (1) define the vertices of a tessellation of the plane 
V a ,b,c with equilateral triangles. 

The above example is saying and it is supporting evidence that the parametrization (2) 
is essentially the only one and the Theorem 1.1 can be extended to cover the degenerate 
solutions case also. 

2 The arguments of the Theorem 1.1. 

Proof. The first part of the Theorem 1.1 follows from [1]. For the second part we are going 
to reconstitute some of the details started in the proof of the particular case: gcd(d, c) = 1. 

Let us start with a triangle in V a ,b,c sa Y AOPQ with P(u ,v ,w Q ) and Q(x ,y , z ). 

By Theorem 4 in [l]we have 

u cv - bw 

Xo = Y 2d 

v aw - cu . . 

y° = ^ ± 2d (6) 

_ w bu - av 
Zo ~Y 2d ' 

for some choice of the signs. This means that d must divide cv — bwo, owq — cuq and 
bu — awo. 

So, we need to look at the the following system of linear equations in m and n: 

Uo = m u m — n u n, 

Vq = m v m — n v n, (7) 
w = m w m — n w n. 

By the Kronecker-Capelli theorem this linear system of equations has a solution if and 
only the rank of the main matrix is the same as the rank of the extended matrix. Since 
m v n w — m w n v = ad, m w n u — m u n w = bd and m u n v — m v n u = cd (one can check these 
calculations based on the definitions in (2)) then the rank of the main matrix is two and the 
rank of the extended matrix is also two because its determinant is u ad + v bd + w$dc = 0. 
This implies that (7) has a unique real (in fact rational) solution in m and n. 

We want to show that this solution is in fact an integer solution. Solving for n from each 
pair of equations in (7) we get 

^ _ v m w - w m v _ w m u - u m w _ u m v - v m u 
ad bd cd 
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Because gcd(a, b, c) = 1, there exist integers a', b' , d such that aa' + bb' + cd = 1. Then 
one can see that 



n _ a'(vom w - w m v ) + b'(w m u - u m w ) + c'(u m v - v m u ) ^ 

Next, from (9) we observe that in order for n to be an integer it is enough to prove that 
d divides v m w — wom v , wom u — uom w and uom v — vom u . Hence we calculate for example 
v Q m w — Wom v in more detail: 

das — rbc v Q qr — (das — rbc)w 

v m w - w m v = v r w = = 

q q 

—daswo + v (3d 2 — c 2 )r + rbcw c(bw — v c)r + 3rv d 2 — dasw 



q q 
From (6) we see that bwo — Vqc = ±d(2xo — uq). Hence, 

d[±c(2x - u )r + 3rv d - asw ] 
v m w - w m v = . (10) 

q 

Assuming that gcd(d, c) = ( we can write d = (d\ and c = (ci with gcd(di, c{) = 1. Also 
we see that ( 2 must divide q = 3d 2 — d 2 so let us write q = ( 2 qi. If p is a prime dividing 
C, it must be an odd prime and if it is of the form Ak + 3 it must divide a and b which is 
contradicting the assumption that gcd(a, b,c) = 1. Therefore it must be a prime of the form 
4k + 1. Hence q\ is still a sum of two squares. 

In the proof of Theorem 13 in [1] one can choose r and s with the extra condition 
that r and s are divisible by (. Indeed, Lemma 14 in [1] is applied to (ad) 2 + 3(db) 2 = 
£ 2 [(aci) 2 + 3(dib) 2 } and to q — ( 2 qi but instead one can apply it to (aci) 2 + 3(dib) 2 and to qi 
giving, let us, say T\ and s\. Then we put r = Qr\ and s = £si and then all the arguments 
there go as stated. From (10) we see that 



v m w - w m v = 

d[±c(2x — u )r + 3rv d — asw ] (di [±(ci (2x — n )C^i + 3Crit>oC^i — aC^i^o] 
q C 2 Qi 

C 2 d 1 [±ci(2x - Mo)^iC + 3riw c/iC - asiw ] _ di£ 

(ii) 

where £ = ±Ci(2x — u )ri(+3riv di(— asiiu . This implies that d\ must divide v m w — w m v . 
In a similar way we can show that d\ divides Wom u — UqVtl w and u m v — v m u and so from 
(9) we see that n is a rational with denominator having £ as a multiple. Similar arguments 
will give us that m is of the same form. 

The triangle having the coordinates as in (1) with these m and n will give the length 



I = 2d (m — ran + n 



2\ 
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But this whole construction can be repeated for a or b instead of c and we obtain that 



I 2 = 2d 2 (m 2 — mini + n 2 ), 

for some rational numbers m 1; ni with denominator having multiple rj = gcd(d, b) for exam- 
ple. Since gcd((, 77) = 1 we see that m 2 — mn + n 2 = m\ — mini + n\ must be an integer. 
Therefore 

I 2 = 2d 2 (a 2 -a(3 + fi 2 )^^ eZ. (12) 
On the other hand if Uq + v 2 + w 2 = I 2 and ( divides d we can see that 

u 2 Q + v 2 + w 2 Q = {mod C 2 ). (13) 

We also know that au + bv + cw = and hence a 2 ul = b 2 VQ + 2bcv w + c 2 w 2 . This 
implies a 2 u^ = b 2 v\ + 2bcv w (mod C 2 )- But a 2 + b 2 = 3d 2 — c 2 = (mod ( 2 ) too and then 
a 2 («o + Vq) = 2bcvoWo (mod ( 2 ) which combined with (13) gives 

a 2 w 2 Q + 2bcv w = (mod ( 2 ). (14) 

Because we must have gcd(a,Q = 1, (14) implies that ( divides w . Indeed, if p is a 
prime that has exponent one in the decomposition of ( then (14) gives in particular w 2 = 
(mod p) and so p must divide w . If the exponent of p in £ is two, then (14) in particular 
implies that wo is divisible by p but then awp = (mod p 3 ) which implies p 2 divides w$. 
Inductively if the exponent of p in ( is k then this must be true for wq too. Hence we must 
have ( a divisor of wo so wq = (w' . 

Now we can go back to (11) and observe that we can rewrite it as 

v r - w Q m v = , (15) 

qi 

where £' = ±ci(2xo — u )ri + 3riVodi — asiiu' . Now we observe that the left hand side of 
(15) is a multiple of ( since r and are - After simplification with ( this will show that 
qi must divide in fact £' and so d divides v m w — w m v . Similar arguments can be used to 
deal with the other cases. Hence n must be an integer and so should be m. Changing the 
variables as in (5), one of the corresponding triangles given by (1) is going to match with 
the triangle OPQ. □ 

Remark: One can see that the condition on r and s to be divisible by ( is implied by 
asking only that the numbers in (2) be integers. Indeed, given such a choice of r and s, they 
will define by (1), in which m = 1 and n = 0, an equilateral triangle with integer coordinates. 
According to the proof of Theorem 1.1, w = r and z = (r — s)/2 must be divisible by (. 
This implies that r and s must be multiples of (. Therefore any parametrization as in the 
Theorem 1.1 is unique up to the transformations (5). 
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Figure 1: The graph of g extrapolating / over the interval [1, 2000] 

3 Behavior of the sequence ET(n) 

The calculations of the ET(n) for all n < 1105 gave us enough data to be able to extrapolate 

the graph of n — » as shown in Figure I. The function we used to extrapolate 

is of the form g(x) = a + having clearly a as limit at infinity. Then we made it 

agree with / on three points.* That gave us a : = 5.079282921, b := -0.7091588389, and 
c := —0.8403164433. Numerically then we discovered that the average of \f(k) — g(k)\ over 
all values of k = 1, 1105 is approximately 0.002638971108. 

One conjecture that we would like to make here is that f(n) is a strictly increasing 
sequence and then as result it is convergent to a constant C ~ 5.08. 

The graph of the "derivative" of ET(n) (Figure 2) is almost like the graph of h(x) = 
C(x + l) k where k := 4.151431798 and C := 2.660972140. The third difference of ET{n) as 
represented in Figure 3 seems to bring a chaotic flavor to this sequence and it is saying in a 
certain sense that no simple formula for ET(n) can exist. 
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Figure 2: The graph of n — > ET(n + 1) — ET{n) 



Third difference sequence of ET(n) 
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Figure 3: A 3 ET(n) 
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(Concerned with sequence A102698 .) 
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